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ABSTRACT 

We combine Spitzer 24fim observations with data from the COMBO-17 survey for ~15,000 0.2 < z < 1 
galaxies to determine how the average star formation rates (SFR) have evolved for galaxy sub-populations 
of different stellar masses. In the determination of <SFR> we consider both the ultraviolet (UV) and the 
infrared (IR) luminosities, and account for the contributions of galaxies that are individually undetected at 
24/zm through image stacking. For all redshifts we find that higher-mass galaxies have substantially lower 
specific SFR, <SFR>/<M*>, than lower-mass ones. However, we find the striking result that the rate of decline 
in cosmic SFR with redshift is nearly the same for massive and low-mass galaxies, i.e. not a strong function 
of stellar mass. This analysis confirms one version of what has been referred to as 'downsizing', namely that 
the epoch of major mass build-up in massive galaxies is substantially earlier than the epoch of mass build-up in 
low-mass galaxies. Yet it shows that star formation activity is not becoming increasingly limited to low-mass 
galaxies towards the present epoch. We argue that this suggests that heating by AGN-powered radio jets is 
not the dominant mechanism responsible for the decline in cosmic SFR since z ~ 1, which is borne out by 
comparison with semi-analytic models that include this effect. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: general 



1. INTRODUCTION 

A key observable statistic that describes the evolution of 
the galaxy population is the average SFR as a function of 
epoch and of galaxy stellar mass. Over the last decade, 
much progress has been made towards delineating, and in 
part understanding, the relationship between star formation 
history and stellar mass. "Archaeological" studies of present- 
day galaxies demonstrate a strong correlation between star 
formation history and present-day stellar mass in the sense 
that the bulk of stars now in massive galaxies must have 
formed at earl ier epochs than stars now in less massive 
galaxies (e.g. , iKauffman net alJ l2003b iThomas et all 12005b 
Galla zzTet alJl2005t iPanter et al.l 12006). Surveys of cosmo- 
logically distant galaxies have revealed a related trend, often 
referred to as "downsizing" in that context. Downsizing most 
commonly refers to the observation that intense star formation 
(M*/SFR <sc t Hubble) becomes increasingly limite d to systems 
of lower and lower mass as redshift decreases dCowie et alJ 
19961: IChaoman et al.l 12001 Uuneau et all 120051: iBauer et all 
20051: iFeulner et al. 2005: lDaddi et alJl2005h iBell et alJl2005h 
Pap ovich et al. 120061) . Recently, this phenomenology has 



been recast in more quantitative terms through study of the 
specific SFR (i.e., SFR per unit stellar mass) as a function of 
redshift, finding that at z < 1, low-mass galaxies appear to be 
forming stars at a higher rate per unit mass than their more 
massive cousins (fBr inchmann et al. 2002 IBauer et all 120051 : 
iNoeske et al1l2007l:lFeulner et al.ll2005l) . 

There are two important parts of this empirical picture that 
require clarification. Firstly, many studies use SFR indi- 
cators which are highly susceptible to d ust extinction (e.g. , 
[On] or UV-optical SED-derived SFRs; ICowie et all 1 19961: 
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Uuneau et all 120051: IBauer et alJ 120051: iFeulner et all 120051) . 
Secondly, all studies to date have relied on individual detec- 
tions to allow SFR derivation. Neglecting SFR contributions 
from galaxies with individually undetectable star formation 
tracers biases the census of SFR at all masses. However, the 
bias is especially severe at low masses where the SFR limiting 
sensitivity implies that all low-mass galaxies with individual 
SFR detections have extremely high specific SFR. 

In this letter, we present an analysis of Spitzer 24 /mi ob- 
servations of two COMBO-17 fields designed to circumvent 
these two key limitations. We estimate SFR by combining 
indicators of obscured and unobscured star formation, and 
attempt to derive average SFRs through stacking of galaxy 
sub-samples with individually-undetected SFR indicators. In 
§2, we summarize the data and galaxy samples. In §3, we 
describe our SFR measurements. In §4, we present the con- 
tribution of different bins in stellar mass to the cosmic SFR, 
and present the evolution of specific SFR with lookback time, 
again as a function of stellar mass. We discuss our results 
in §5. We adopt a cosmology with Ho =70kms _1 Mpc~', 
H M = 0.3 andO A = 0.7. 

2. THE DATA AND GALAXY SAMPLES 

The COMBO-17 survey dWolf etalj 120031) has imaged 
three 30' x 30' fields in five broad and 12 intermediate op- 
tical bands from 350 to 930 nm with the goal of deriving 
high-quality photometric redshifts 4 for > 1 000 galaxies with 
ntR < 24 mag in each field. The rest-frame luminosities in the 
U, B and V band, and a synthetic UV2800 band centered at 
280 nm have been derived for all COMBO-17 galaxies with 
photo-z's. By fitting the COM BO-17 multi-band p hotometry 
with a library of SED templates. iBorch et alJ(l2006l) estimated 
the galaxies' stellar mass-to-light ratio (M/L), a nd hence stel- 
lar ma ss for all otr < 24 galaxies based on a Krou paet alJ 
d 1993b initial mass function (IMF; such an IMF has stellar 



4 Redshift unc ertainties are fe/( l + z) ■ 
tude of m R ~ 22 ( Wolf et all 20041) . 
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masse s similar to wi thi n 10% compa red to those derived us- 
ing alChabrierll2003l or lKroupall200ll IMF). From these data, 
iBorch et all (120061) constructed stellar mass functions in four 
even redshift slices between z = 0.2 and z — 1 . 

Two of three COMBO- 17 fields, the extended Chandra 
Deep field South (CDFS) and the Abell 901/902 (A901) 
have been observed by the Multiband Imaging Photome- 
ter on Spitzer (MIPS: iRieke et al.ll2004l) at 24yt/m. The re- 
duced 24 yum images cover a n area of 1° x 0°5 around the 
CDFS (Pa povichetaTl l2004T) and of 30' x 55' in the A901 
field dBell et al.ll2007l) . Both images have a pixel scale of 
l'/25 pixel -1 , a point spread function (PSF) with full width 
at half maximum (FWHM) of ^6". The IRAF/DAOPHOT 
package and empirically-determined PSFs are used to simul- 
taneously fit multiple sources to the 24 yum images and derive 
total flux for each detected object to a depth of 83 //Jy at the 
5 cr le vel (to a completeness of 80%; see also iPapovich et all 
2004 for details of data reduction, source detection and pho- 
tometry). To match these 24/mi catalogs to the COMBO-17 
galaxy samples, a position tolerance of 2'.'2 is adopted. If 
more than one galaxy exist within the tolerance, the near- 
est one is selected. The overlap area between COMBO- 
17 and MIPS imaging contains 9785 and 12995 galaxies at 
m R < 24 in the CDFS (~800 square arcminutes) and the A901 
field (~850 square arcminutes), respectively. Of them, 1735 
(1872) are individually detected at 24yum. Only <4% (6%) of 
the 24 yum-detected sources, namely 64 (108), have multiple 
(mostly two) COMBO-17 galaxies within the tolerance. We 
select galaxies with stellar mass M, > 10 y M in the redshift 
range 0.2 < z < 1, resulting in a sample of 6893 (8184) galax- 
ies in the CDFS (A901) field. Of the 6893 (8184) galaxies, 
1201 (1250) are individually detected at 24yum. The average 
24 /im flux of the galaxies individually undetected by MIPS 
will be estimated by stacking, as described in the next section. 
The samples are complete at ma ss limit ~1, 1.8, 3, 6xlO I() M 
for redshift z -0.3, 0.5, 0.7, 0.9 (IBorch et al.ll2006l) . 

Owing to mismatches in depth between wide-area IR sur- 
veys (i.e., IRAS) and well-characterized optical galaxy sur- 
veys (i.e., the SDSS and 2dFGRS), it is challenging to pro- 
duce an appropriate and well-defined local comparison sam- 
ple. We have attempted to produce an adequate local con- 
trol sample by collecting a sample of 2177 galaxies from 
the NASA/IPAC Extragalactic Database (NED) with 2MASS 
TT-band magnitude K < 12 in the volume of 1500 < 
cz/(kms -1 ) < 3000 and Galactic latitude b > 30°. Assum- 
ing a /f-band M/L of O.6M O /L and a Kroupa IMF, stellar 
masses were derived from the /^-band absolute magnitude. 
The SFR is estimated from the total IR luminosity, as deter- 
mined from IRAS observations. The local sample is not a 
complete volume-limited sample. However, we estimate that 
above M, > 10 10 M o the sample is ~65% complete with an 
incompleteness which is largely geometric in origin (because 
some areas of the sky were not spectroscopically surveyed as 
thoroughly as others). If indeed the incompleteness is primar- 
ily geometric in origin, the avera ge SFRs at a g iven stellar 
mass should not be biased. 5 See iBell et all d2005l) for more 
details about the local sample selection and completeness. 

3. SFR DETERMINATION 

We split the sample galaxies into four redshift slices evenly 
covering the range from z = 0.2 to z = 1. In each red- 

5 At any rate, inspection of Fig.[T]and Fig.f2]shows that none of our con- 
clusions actually depend on the z, ~ data; the purpose of its inclusion is 
simply to extend the trends to z ~ 0. 



TABLE 1 

The cosmic star formation history split by galaxy stellar mass. 
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shift slice, sample galaxies are divided into five stellar mass 
bins: four bins with a bin width of 0.5 dex spanning from 
M« = 10 9 M G to 10 n M G , and one bin for M« > 10 n M o . 
Each sub-sample of galaxies defined by optical data contains 
individually-detecte d and indiv i dually -undetected objects at 
24 /im. Following Zhen g~et~aT] j2006h we estimated the to- 
tal 24 yum flux for the sub-sample of individually-undetected 
objects through stacking. The total 24ytmi flux of the sub- 
sample is then obtained by adding the individual fluxes from 
the detected objects and stacked flux from the individually- 
undetected objects. Our SFR indicator, the total IR luminosity 
(8-1000 yum), is extrapolated from the 24 yum luminosity us- 
ing thr e e sets of local IR SEP tem plates from | Cha rv & Elbaz 
(120011) iDale & Heloul (120021) and lLagache et alJ (120041) . In 
IZheng et all (120071) . we use stacking of IR-bright galaxies 
with 0.6 < z < 0.8 at 70/im and 160 /mi to show that the 
average IR SED of an IR-bright galaxy is spanned by these 
local template sets, lending credibility to the total IR lumi- 
nosity estimated here to within a factor of two. 

We use the bolometric luminosity produced by young stars 
(UV+IR) to estimate the SFR, a procedure that accounts for 
both obscured and unobscured star formation. The total UV 
luminosity (1216A-3000 A) is estimated from the COMBO- 
17 2800A rest-frame luminosi ty L v ,2800 as Lpy = 1.5 vL Vj 2800- 
We used the conversion from Bell et al. ( 200 5b to convert the 
UV+IR into SFR. The SFR is calibrated to a lChabrierl d2003l) 
IMF. Using the sum of stellar masses for each sub-sample 
from IBorch etail ((2006), we calculated the average specific 
SFR, i.e. < SFR > / < M« >, for all redshift and stellar 
mass bins. Formal uncertainties were estimated from boot- 
strapping. For each redshift and stellar mass bin, the average 
between the CDFS and A901 was taken as the final result; the 
difference between the two is taken as a measure of field-to- 
field variance and included into the uncertainties. The esti- 
mates of both SFRs and stellar masses are affected by system- 
atic errors, that other studies have estimated to be < 0.3 dex. 
These systematic uncertainties are not included in our error 
bars, but are substantially smaller than the dynamic range of 
the trends discussed in this Letter. Due to the COMBO-17 R- 
band selection (tor < 24), our sample galaxies are incomplete 
for low-mass bins at z > 0.4, where low-mass red galaxies 
with little blue radiation from newly-formed stars fail to make 
the /?-band selection cut but low-mass blue galaxies contain- 
ing more on-going star formation remain available for the se- 
lection. The selection effect may lead to a potential overesti- 
mate of the specific SFR for these incomplete bins. The spe- 
cific SFR at z ~ was estimated using the local sample only 
for three high-mass bins where the sample is representative. 

4. <SFR> (Z) AS A FUNCTION OF GALAXY STELLAR MASS 
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Fig. 1 . — The volume-averaged SFR as functions of redshift, for galaxy sam- 
ples of differing mean stellar masses. The average SFR is derived from the 
average bolometric luminosity (UV+IR) for galaxy subsets limited in mass 
and redshift. The gray symbols are existi ng measurement s of the overall 
SFR density in the literature, adapted from Hopkins 12004). The solid cir- 
cles show the overall SFR density derived from the SFR functions given in 
IBell et alj 420071) . The data imply that the fraction of the total SFR con- 
tributed by galaxy sub-samples of different stellar masses has not changed 
since z ~ 1, despite the factor ~7 decline in the total SFR in that interval. 



In the previous section, we have determined the average 
SFRs for galaxies split into stellar mass and redshift bins. In 
order to determine the contribution of different stellar mass 
bins to the cosmic (average) SFR, we must combine the aver- 
age SFRs at a given mass with the space density of galaxies 
in that same mass range, as a function of redshi ft. For this 
purpose, we adopt the stellar mass funct i ons of iBorch etafl 
(2006) for 0.2 < z < 1, and IBell et ail d2003l) for z ~ 0, 
although the results are unaffected to within th e uncertain- 
ties if we adopt stellar mass fun ctions from, e . g..lbrorv et ail 
(12004 . iFontana et all (120061) or lBundy etaT] (120061) instead. 
We calculate the volume-averaged SFR over the redshift range 
< z < 1 for three galaxy stellar mass bins: M»/M > 10 , 
10 10 - 10 11 and 10 9 - 10 ll) . The stellar mass functions from 
IBorch et all d2006l) are derived from -25 000 COMBO-17 
galaxies selected from three separated 30' x 30' fields. There- 
fore, the field-to-field variance is somewhat suppressed in this 
way (as one can see from compariso n of the total SF R, de- 
rived from two COMBO-17 fields bv lBelletal] d2007l) com- 
pared to the smooth trends followed by our mass-binned cos- 
mic SFRs). The local data point is not available for the stellar 
mass bin 10 9 < M*/M < 10 10 . 

The results are listed in Table Q] and shown in Fig.Q] along 
with current measure ments of the co s mic SFR density avail- 
able in the literature dHopkinsll20"04 . IBell et alJ (120071) con- 
structed SFR functions using the same data sets and SFR esti- 
mator (UV+IR) as adopted here. The data points given by 
their SFR functions are shown also in this plot. Our esti- 
mate of the cosmic SFR density, i.e. the sum of the volume- 
averaged SFRs from the three stellar mass bins (contribution 
from galaxies of stellar mass < 1O 9 M is negligible), is in 
good agreement with Bell et al.'s and other measurements 
from the literature. 

Fig. Q] shows that the volume-averaged SFR of galaxies 
in each of the three stellar mass bins decreases towards the 
present day with an indistinguishably similar slope since z = 
1 . In other words, the contributions to the cosmic SFR den- 
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Fig. 2.— The specific SFR (SFR per unit stellar mass), < S FR > / < M, >, 
as function of redshift for galaxy samples of differing stellar masses. The 
average SFR is derived from the average bolometric luminosity (UV+IR) for 
galaxy subsets limited in mass and redshift. The data points in the hatched 
region are incomplete due to the COMBO-17 source selection (/br < 24). 
The gray strip illuminates the overall specific SFR, i.e. the ratio of the cosmic 
SF R density to the co smic stellar mass density, following |Hopkins| 420041) 
and Borch et al. (2006). Local data points are estimated for three high-mass 
bins using a local galaxy sample with near-IR and IRAS observations. The 
errorbar represents the typical uncertainty of our specific SFR determinations. 
The dark-gray dashed line is the specific SFR required to double the mass of a 
galaxy by the present day. Data points above the mass-doubling line represent 
galaxies in intense star formation mode and those below the line represent 
galaxies in quiescent star formation mode. For all galaxy sub-samples of 
different stellar masses, the specific SFR increases with redshift out to z ~ 1 ■ 
At all redshifts, more massive galaxies always have lower specific SFR than 
less massive galaxies. These relations delineate a picture in which lower-mass 
galaxies transition from an intense into a quiescent star formation phase since 
Z ~ 1, whereas most massive galaxies have formed stars quiescently across 
this epoch. Open symbols show results from Feulner et al. 12005) based on 
SFR indicator rest-frame UV luminosity for stellar mass bins (converted to 
a Kroupa IMF) 10 8 25 < M,/M Q < W 925 (triangles), W 925 < MJM e < 
10 10 - 25 (squares) and 10 1025 < M,/M < 10 1L25 (circles). 

sity from galaxies of different stellar mass ranges has not 
changed with cosmic epoch during the last 8 Gyr. Galaxies 
in the 10 10 < M,/M < 10 11 contain the bulk of the overall 
star formation at all cosmic epochs which we have examined. 
These galaxies also dominate the cosmic stellar mass density 
(IBorch et alj|2006l) . A smaller part of the overall star forma- 
tion takes place in galaxies with 10 9 < M,/M < 10 10 . Rela- 
tively little star formation is associated with massive galaxies 
withM. > 10 n M o . 

Fig- Elpresents our measurements of specific SFR, i.e. SFR 
per unit stellar mass. Data points of the same stellar mass 
bin are linked together, showing the redshift evolution of the 
specific SFR. Data points of the same redshift delineate the 
specific SFR as a function of stellar mass. Those enclosed 
in the shaded region are in bins affected by serious incom- 
pleteness due to the COMBO-17 flux limit. For compari- 
son, we overplotted the average specific SFR for all galax- 
ies, i.e. the ratio of the cosmic SFR density to the c osmic 
stellar mass densi ty, calculated using Hopkins (2004 and 
IBorch etafl (120061) . As shown in Fig. [2] the average spe- 
cific SFR for a galaxy stellar mass bin increases rapidly with 
redshift, in qualitative agreement with previous studies (e.g., 
iBauer et al. 20051 IFeulner et ai1l2005h . Compared to our re- 
sults, SFR estimates from rest-frame U V luminosity (open 
symbols in Fig. |2j IFeulner et al.l 12005) produce lower spe- 
cific SFRs for high-mass galaxies and higher specific SFRs 
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for low-mass galaxies due to the reasons mentioned in §Q] 

Yet, the change in specific SFR is nearly equally strong for 
all galaxy mass bins. At all redshifts, the population-averaged 
specific SFR of less massive galaxies is higher than that of 
more massive galaxies. Note that, as expected, galaxies of 
stellar mass around 10 < M,/M Q < 10 11 have specific 
SFR comparable to the 'cosmic average' specific SFR, as they 
contribute most to both the cosmic average SFR and stellar 
mass budget. 

5. DISCUSSION AND SUMMARY 

We combine Spitzer 24 pm observations with COMBO-17 
data to explore the evidence for the 'downsizing' of the star 
formation activity since z ~ 1 . Specifically, we estimate SFR 
from the UV (for unobscured star formation) and the IR lumi- 
nosities (for obscured star formation) for a sample of ~15 000 
galaxies selected from the CDFS and A901 fields. Through 
image stacking, we account for the SFR contributions of oth- 
erwise undetectable galaxies. Our measurements of average 
SFR for galaxy subsets limited in mass and redshift account 
for contributions from both star-forming galaxies and quies- 
cent galaxies in each subset. We find that the relative drop 
in SFR from z ~ 1 to the present day is nearly identical 
for different sub-samples that span a range of ~100 in stel- 
lar mass. That is, all galaxy samples, independent of mass, 
exhibit the same gradual decline of star formation (and spe- 
cific star-formation) over the last 8 Gyr. At all redshifts, we 
find that high-mass galaxies have substantially lower specific 
SFR than their low-mass cousins. 

While many of the qualitative feature s of this paper have 
been understood for some time (e.g.. princhmann & Ellis 
l2000HBauer et al.ll2005l:lJuneau et al.ll2005b . this Letter quan- 
tifies the relationships between specific SFR, galaxy stellar 
mass and redshift. Archaeological' downsizing — that the 
stars in massive galaxies formed at earlier times than the stars 
in low-mass galaxies — is clearly recovered in this paper, in 
Fig. [2] On average, low-mass galaxies have specific SFRs ca- 
pable of more than doubling their stellar mass over the Hubble 
time (at all redshifts), whereas galaxies with M, > 10 11 M 
are not capable of such growth at z < 1 from in situ star for- 
mation alone. 

Arguably the most striking and new contribution of this 
Letter is that the decline in SFR from z ~ 1 to the present 
occurs at the same rate in galaxies of all masses. A num- 
ber of mechanisms have been postulated to drive the demise 
of star formation since z ~ 1: AGN feedback, the declining 
galaxy merger rate, gas consumption, changes in the avail- 
ability of cooling or infalling fresh gas, and a host of other 
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possibilities. We speculate that the declining supply of new 
gas falling in from the cosmic web at z < 1 seems an attrac- 
tive candida te for driving much of the phenomenology that 
we observe dNoeske et al.ll2.007h . The roughly constant rate 
of decline of SFR and SSFR for galaxies of widely differing 
masses is in contrast to the predictions of theoretical mod- 
els in which cooling and hence star for mation is quenched 
by heating due to radio-loud AG N (e.g.. ICroton et al.ll2006t 
iBower et all 120061: but see also iNeistein et al .1120061) . This 
process is more important in more massive galaxies for two 
reasons. First, more massive galaxies are more likely to host 
massive black holes, and more massive black holes are gen- 
erally assumed to be capable of producing more energy and 
stronger heating. Second, it is generally assumed in such 
models that AGN can only heat gas that is cooling from a 
quasi-hydrostatic hot halo, which form in large mass halos 
(MhMo > ~10 I2 M o ) as opposed to the "cold mode" flows that 
ar e dominant in lowe r mass halos. For example, the models 
of lBower et al.1 (120061) predict a decline of a factor of ~8 from 
Z ~ 1 to the present, in the SFR contributed by galaxies more 
massive than 1O 1O 8 M (see their Figure 8), while the SFR 
contributed by less massive galaxies declines by a factor of 
about 2-4 over this same interval. We find, using an indepen- 
dent model with 'radio-mode' heating by AGN (Somerville 
et al. in preparation), that the specific star formation rate in 
galaxies with stellar mass greater than 10 11 M has declined 
by a factor of ~15 since z ~ 1, while the specific SFR in 
the lower mass bins (the same bins used here), declines by 
only a factor of 3-5. The precise rate of decline of the SFR 
as a function of galaxy mass is no doubt sensitive to the de- 
tails of the physics of star formation, but the prediction of a 
strongly differential 'quenching' of star formation by AGN 
feedback in the form of radio heating, as it is currently imple- 
mented, seems inevitable. Our results therefore suggest that 
radio mode heating by AGN is not the dominant mechanism 
responsible for the decrease in cosmic SFR from z ~ 1 to the 
present. 
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